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ABSTRACT 

As a logical next step in improving the radio-to-submm spectral index as a 
redshift indicator (Carilli & Yun), we have investigated a technique of using the 
entire radio-to-FIR spectral energy distribution (SED) for deriving photometric 
redshifts for dusty starburst galaxies at high redshift. A dusty starburst SED 
template is developed from theoretical understanding on various emission mech¬ 
anisms related to massive star formation process, and the template parameters 
are selected by examining the observed properties of 23 IR selected starburst 
galaxies: = 58 K, /? = 1.35, and fnth = 1- The major improvement in using 

this template SED for deriving photometric redshifts is the signihcant reduction 
in redshift uncertainty over the spectral index technique, particularly at higher 
redshifts. Intrinsic dispersion in the radio and FIR SEDs as well as absolute 
calibration and measurement errors contribute to the overall uncertainty of the 
technique. The derived photometric redshifts for Eve submm galaxies with known 
redshifts agree well with their spectroscopic redshifts within the estimated un¬ 
certainty. Photometric redshifts for seven submm galaxies without known spec¬ 
troscopic redshifts (HDF850.1, CUDSS14.1, Lockman850.1, SMM J00266-I-1708, 
SMM J09429+4658, SMM J14009+0252, FIRBACK J1608+5418) are derived. 

Subject headings: galaxies: high-redshift - galaxies: starburst - infrared: galax¬ 
ies - submillimeter - radio continuum: galaxies - techniques: photometric 


1. Introduction 


Sensitive observations at snbmm wavelengths are revealing what may be a population 
of active star forming galaxies at high redshift which are unseen in deep optical surveys due 
to dust obscuration (Small et al. 1997; Barger et ah 1998; Hughes et ah 1998; Bales et al. 
1999). The differential source counts clearly indicate a large excess of far-infrared sources 
by a factor of 10-50 over the no-evolution models derived from the optical deep survey data 
(Guiderdoni et al. 1998; Blain et al. 1999; Matsuhara et al. 2000; Scott et al. 2001), and a 
large fraction of star formation may be hidden by dust. The analysis of the snbmm SCUBA 
source counts and inferred redshift distribution based on the radio-submm flux density ratio 
(Carilli & Yun 1999) suggests that luminous dusty galaxies may dominate the star formation 
history at early epochs (;2 ~ 1 — 3) (Barger et al. 2000). Many of these faint snbmm sources 
are found to be very faint, red {R >25,K> 21) sources (Small et al. 1999), whose optical 
redshifts may be inaccessible even for the 10-m class telescopes. These galaxies are nearly 
entirely missing from the optical studies of star formation at high redshifts (e.g. Steidel et 
al. 1999; Chapman et al. 2000), and a signihcant revision to the optically derived cosmic star 
formation history may be needed (see Blain et al. 1999). 

Because direct redshift measurements are not possible in most cases, the redshift dis¬ 
tribution and the cosmic evolution of the dusty snbmm galaxy population are not well 
determined at the moment. In a recent paper, we have proposed a technique of using the 
radio-to-submm spectral index as a redshift indicator. This technique is based on the uni¬ 
versal radio-to-far infrared (FIR) correlation for star forming galaxies (Condon 1992), with 
the assumption that the spectral shapes may be similar enough to be able to differentiate 
between low and high redshift objects (see Carilli & Yun 1999). To understand the magni¬ 
tude of scatter in this relation and resulting uncertainty in the redshift estimates, we derived 
an empirical af™ — 2 : relation^ using the observed spectral energy distributions (SEDs) of 
17 low redshift starburst galaxies (Carilli & Yun 2000; also see Dunne, Clements, & Bales 
2000a). There is a signihcant scatter, but the existing data on high redshift star forming 
galaxies and AGN-hosts appear to follow this relation well. 

While the radio-to-submm spectral index has been demonstrated to be a useful redshift 
indicator, using it to derive a redshift for any particular object may be risky because of the 
scatter in the observed q;J ^4 — relation and the battening in the relation at high redshifts. 
Obtaining redshift estimates using hux density ratios at other wavelength bands has also 
been suggested although none is particularly more successful (e.g. Hughes et al. 1998; Blain 
1999; Hines & Low 1999; Fox et al. 2001). More accurate redshift estimates for dusty 
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starbursts may be obtained by utilizing the information associated with the entire radio-to- 
FIR SED. Such a photometric redshift technique requires a good SED template, and here 
we derive an SED template based on the theoretical expectations of thermal dust, thermal 
Bremsstrahlung, and non-thermal synchrotron emission from a dusty starburst galaxy. Our 
SED template is constructed by examining the properties of the 23 IR selected starburst 
galaxies whose submm/EIR SED data are readily available in the literature. When applied 
to a sample of submm galaxies with known redshifts, the resulting redshift estimates are in 
excellent agreement with the spectroscopic redshifts. We also report photometric redshifts for 
a sample of submm galaxies without known spectroscopic redshifts, including the brightest 
SCUBA source in the Hubble Deep Field (Hughes et ah 1998). 


2. Starburst SED Model 

The observed SED from radio to infrared wavelengths for a dusty starburst galaxy is 
dominated by the energetic signatures of young stars and their environment. The bulk of the 
radiation from young OB stars is absorbed by dust and re-emitted in the infrared. The HII 
regions surrounding the young stars are also the sources of thermal Bremsstrahlung (free- 
free) emission. Energetic electrons produced in supernovae and their remnants are thought to 
account for nearly all of the radio synchrotron radiation. Following the discussion by Condon 
(1992), a simple starburst SED model is constructed as a linear sum of dust continuum {Sd), 
thermal Bremsstrahlung or free-free (S'th), and non-thermal synchrotron {Snth) emission. 


2.1. Thermal Dust Emission 

The flux density, Sd{n), from a cloud at a distance containing N spherical dust 
grains each of cross section temperature T^, and emissivity Q{n), is given by 

Sd{u)=N{ad/Dl)Q{u)B{u,Td), (1) 

where B{u,Td) is the Planck function describing thermal black body radiation. The term 
N{ad/D\) is the solid angle VLd subtended by the dust source in the sky. For the emissivity 
function, we adopt a form 

Q{v) = I - exp [-{ — )% (2) 

where (3 is the dust emissivity index generally thought to be between 0 and 2 (Hildebrand 
1983). This particular functional form allows the dust spectrum to be that of a purely black 
body above the critical frequency Vc where the dust clouds become optically thick while it is 
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that of the grey body spectrum at lower frequencies. Observations suggest that thermal 
dust emission is nearly optically thick even at A ~ 100 pm for galaxies dominated by intense 
starburst regions (see Scoville et ah 1991; Solomon et ah 1997), and we adopt Uc = 2 x 10^^ 
Hz (Ac = 150 pm). Combining the Eqs. 1 & 2, 

5,(z/) = 0,H(z/,T,)(l-e-(^)"). (3) 

For a dust source of 9 arcsec in diameter, expected flux density at z/ (in GHz) due to thermal 
dust emission is 

S,iu) = 2.8 X 10~^ ^;o.o48./r,]_i (l - (4) 

The FIR luminosity, Lpjji, is derived by integrating Fq. 4 over the emitting area Qd and 
over the frequency range corresponding to A = 40 — 500 pm, Lpm = 47rZl| f f Sd{v)di'dVL. 
The massive star formation rate (SFR) for a dusty starburst galaxy is then directly re¬ 
lated to its dust continuum spectrum by the FIR luminosity (i.e. SFR = Lpip/{5.8 x 

lO^L©) Mq yr“\ Kennicutt 1998) 


2.2. Thermal Bremsstrahlung (free-free) Emission 

Thermal Bremsstrahlung (free-free) emission from ionized gas in the HII regions sur¬ 
rounding hot yonng stars can be described in the same functional form as the thermal dust 
emission as in Eq. 3, 

S„{X = - e-^o), (5) 

where Tg is electron temperature and Tff is free-free optical depth at the observed frequency. 
In the Rayleigh-Jeans regime {hu kT), B{u,Tf,) = 2A;TeZ/^/c^ at radio freqnencies, and 

Sff{u) = 2kTyc-^Qff{l-e--ff). ( 6 ) 

Free-free optical depth is given by Tff = 0.083EMu~‘^'^T~^'^^, where EM = J n^dl 
(cm“® pc) is the emission measnre. In the optically thin regime, thermal Bremsstrahlung 
emission then has the familiar dependence. A circular source of 9 arcsec in diameter 
would have flux density of 


Sff{u) = 4.8 X lO-®EMz/-°'iTg-°'3502 jy_ 


( 7 ) 
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2.3. Non-thermal Synchrotron Emission 


For an ensemble of relativistic electrons with a power-law distribution of energy N{E) = 
NoE~'^ and an isotropic velocity distribution, their synchrotron emission coefficient is 

e oc , (8) 


where 6 is the angle between the magnetic held B and the line of sight to the observer. The 
resulting power-law spectrum has a form Si^ oc where a = (7 — 1)/2 is the spectral index. 
Bell (1978) proposed the acceleration of cosmic rays by supersonic shocks (of velocity Vg) 
as the main mechanism for non-thermal synchrotron emission from supernova remnants and 
derived the volume emissivity as, 


e(z/) ~ 3 X 10 


— 33^—0 


n 


a 


( 7 —i)(77F7)( 


\4o/ 


B 


cm~ 




erg s“^Hz-^cm-^ 


(9) 


Assuming supernova remnants dominate the non-thermal synchrotron emission from 
galaxies, Condon (1992) derived a relation Snth i'sn where usn is the Type II supernova 
rate per year (see his Eq. 17). However, this relation predicts a vsn more than an order 
of magnitude too large for our Galaxy, and Condon proposed a revised relation based on 
Galactic normalization (i.e. his Eq. 18), which can be re-written as 

Snthi.^) = 1150 p~'^usnDl~‘^ Jy, ( 10 ) 

where Dl is luminosity distance in Mpc. 


2.4. Starburst SED Template 

The hrst and the foremost important step in all photometric redshift technique is es¬ 
tablishing a robust and reliable template. Key SED template parameters for the starburst 
SED template are and (3 for the dust emission (Eq. 4), Te and EM for thermal free-free 
emission (Eq. 7), and vsn and a for the non-thermal synchrotron emission (Eq. 10). These 
parameters and their dispersions are determined by deriving best £t SED models for the 23 
IR selected starburst galaxies whose submm/FIR data are readily available in the literature. 

Flux density for dust emission Sd can be derived directly using Eq. 3 by specifying the 
source size dust temperature Td, and emissivity (3. The effective source solid angle Qd 
for a given source is derived from the observed 100 pm flux density, assuming dust emission 
is nearly optically thick (see § 2.1). For example, the prototypical ultraluminous galaxy 
Arp 220 is a 114 Jy (1.14 x 10“^^ erg cm“^ Hz“^) source at 100 pm, which translates to 
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= 3 X 10“^^ rad‘^ for = 59 K (see Table 1). This corresponds to a circular area with 1.2" 
(390 pc) in diameter^, which agrees well with the size of the molecular gas complex fueling 
the nuclear starburst traced in CO (Scoville, Yun, & Bryant 1997; Downes & Solomon 1998; 
Sakamoto et ah 1999). The best £t dust SED model is then selected by examining the 
parameter space for Td and (3, between 25-85 K and 1.0-2.0, respectively. 


Thermal Bremsstrahlung or free-free flux densities can be derived from Eq. 7 by de¬ 
termining the frequency where free-free opacity r/j becomes unity and estimating the size 
of the emitting region. Alternatively, one can also compute the free-free flux density from 
the inferred SER as it is proportional to the production rate of Lyman continuum photons. 
Using the Ha normalization of Kennicutt (1998) for a Salpeter IMF and assuming 50% of 
Lyman continuum photons are quenched by dust absorption, Eq. 23 of Condon (1992) can 
be re-written as 

Lffiu) = 7.9 X 1026zz-0-i(^^^^) erg s-^Hz-h (11) 

The free-free flux density for a galaxy with a given FIR derived star formation rate is then 

Q p f? 

Jy. (12) 

The free-free flux density Sff determined this way is nearly identical to the value derived 
using Eq. 7 assuming r// ~ 1 at z/ = 1 GHz and D// = fid- 

Non-thermal synchrotron flux density Snth can also be parameterized as a function of 
SER since usn oc SFR. Adjusting Eq. 20 of Condon (1992) for the Salpeter IMF, Eq. 10 
can be re-written as. 


Snthiy^ — 25 fnth^ 


SFR 


'M^yr 


Jy. 


(13) 


The synchrotron spectral index a is known to he within a narrow range around 0.7-0.8, and 
we simply adopt a = 0.75 to minimize the number of free parameters in our model. Condon’s 
original derivation arbitrarily employed a Galactic normalization as noted earlier. We thus 
add a scaling factor fnth (of order unity) here in order to determine the normalization more 
suitable for starburst galaxies. 


Only the SED data at frequencies near and below the dust peak are used for the model 
fit because the mid- to near-IR data points clearly require additional, higher temperature 
dust components. The SED fit for Arp 220 is shown in Figure 1, and it demonstrates the 
effects of changing Td while holding (3 constant. The scatter of data points about various 


^An ensemble of distributed sources with the equivalent total area is also allowed by this simple model. 
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models suggests that there are systematic scaling differences among different flux density 
measurements, generally larger than the nominal uncertainties reported, especially in the 
submillimeter wavelengths. At a glance, all three dust temperature models plotted seem 
to do an adequate job of htting the submm and FIR data points qualitatively. However, a 
closer examination reveals that the = 74 K model fails to match the submm points. The 
Td = 43 K model matches the submm data points reasonably well, but it clearly falls short 
on predicting the FIR data points. Previous studies of submm dust properties for luminous 
infrared galaxies generally favored cold dust temperature between 40 and 50 K for Arp 220 
(e.g. Bales et ah 1989; Scoville et al. 1991; Dunne et al. 2000b). If we were to choose a single 
temperature dust SED model that best hts all observed data points between 1mm and 60 
pm, however, a slightly warmer dust model is favored. 

Our starburst SED model analysis is applied to 23 IR-selected starburst galaxies that 
are selected for (1) their Lpm exceeding 10^^ L©; and (2) having at least two measurements 
covering their submm part of the spectrum. The SED data are primarily drawn from the 
IRAS Faint Source Catalog, NRAO/VLA Sky Survey (Condon et al. 1998), Condon et 
al. (1991), Carico et al. (1992), Rigopoulou et al. (1996), Benford (1999), Lisenfeld et al. 
(2000), Dunne et al. (2000b), and Dunne & Bales (2001). As summarized in Table 1, the 
characteristic dust temperature ranges between 46 K and 74 K, with a mean of 58 ± 9 K. 
Dust emissivity (3 also ranges widely, between 1.05 and 1.70, with a mean of 1.32 ±0.17 and 
a median of 1.35 (see Figure 2). It is likely that all galaxies consist of an ensemble of dust 
clouds with a range of temperature (see Yun & Scoville 1998; Frayer et al. 1999; Dunne & 
Bales 2001), and single temperature SED models tend to favor a smaller (3 (flatter submm 
dust spectrum) and higher dust temperature than typically found in giant molecular clouds 
(10-20 K). Since available SED measurements are generally sparse in practice, we accept 
these limitations of a single temperature dust model for the purpose of keeping the number 
of free parameters manageably small. 

Thermal Bremsstrahlung (free-free) emission makes signihcant contribution only in the 
bottom of the SED trough between the non-thermal synchrotron and thermal dust feature, 
and it plays essentially no role in dehning the starburst SED template. Variations in the non- 
thermal synchrotron emission are tracked by a parameter f^th (see Eq. 13). The dispersion 
in the radio-FIR correlation is measured to be about 0.25 in a logarithmic scale (e.g. Yun 
et al. 2001), and a rather broad range of fnth found (see Figure 3) is consistent with this 
expectation. Presence of a radio AGN (e.g. Mrk 231, NGC 6240) can account for the increase 
in some cases, but synchrotron emission efhciency may indeed vary from galaxy to galaxy. 
The median for the whole sample is 1.1, but it is reduced to 1.0 if the two clear radio AGN 
hosts are removed. Therefore, on average, the non-thermal synchrotron emission from these 
dusty starbursts appear to follow the Galactic normalization adopted by Gondon (1992). 



In summary, we adopt = 58 K, /? = 1.35, and fnth = 1-0 for our photometric redshift 
SED template. This relatively flat (3 is in good agreement with the analysis of 102 galaxies by 
Dunne et ah (2000b), but the mean dust temperature we adopt is signihcantly warmer than 
the mean of = 36±5 K Dunne et ah derived using 60 p, 100 /i, and 850 /i measurements. 
The main difference is that the dust emission from our IR selected, luminous dusty starbursts 
are on average signihcantly warmer than those of the Dunne et al. sample, which includes a 
large number of late type held galaxies. 


3. Photometric Redshifts for Dusty Galaxies at High Redshift 

3.1. Procedure 


Once the dusty starburst SED template is chosen, then the entire radio-to-FIR contin¬ 
uum spectrum for any given galaxy can be described in terms of just the total star formation 
rate SFR and the luminosity distance Dl. The dust spectrum Sd in Eq. 4 can be re-written 
as 


Rrf(z/) = 1.5 X 10 


_g SFR{l + z)u^ 


D'U 


30 . 00083 !/ 


1 ) 


(1 


-( ^ 

3 '- 2000 / 


Jy- 


(14) 


Combined with the expressions for free-free and non-thermal synchrotron emission given in 
Eqs. 12 & 13, the entire radio-to-FIR continuum spectrum can be written as. 


125 


where SFR is in Mq yr“\ Dl in Mpc^, and z/q = (1 + z)iyobs in GHz. A (1 -|- 2 :) term 
is needed for sources at cosmological distances in order to account for the frequency or 
wavelength folding by the Doppler ehect. 

The redshift and the FIR luminosity {SFR) of a particular galaxy can be determined 
simultaneously by comparing the observed SED with the dusty starburst template in the 
{z,SFR) space. To be precise, one needs to include a correction term accounting for the 
cosmic microwave background (CMB) because most SED data are measured in contrast to 
the CMB. In practice this correction can be safely ignored as long as Td/(1 + ;2) ^ 2.7 K. 

As a galaxy is placed further and further away, its entire SED shifts to the bottom 
(fainter, because of and to the left (lower frequency, Pobs = ^'o/(l + z)) in the {z,SFR) 


^We adopt Ho = 75 km sec ^ Mpc ^ and Qo = 0.5 for this paper. 
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space. This generic behavior for a dusty starburst and the resulting change in the apparent 
spectral index between 1.4 GHz and 850 /iui has been pointed out previously as a redshift 
indicator by Carilli & Yun (1999). The slope of the rising part of the dust spectrum is such 
that the Doppler shift of the spectrum nearly offsets the drop in flux density, making 
the submm bands particularly attractive for blind searches, but the SED measurements on 
this part of the dust spectrum offer limited redshift information for the same reason (e.g. 
Hughes et al. 1998; Fox et ah 2001). As in other photometric redshift techniques, the redshift 
information comes from distinct spectral features such as the trough between the declining 
non-thermal synchrotron emission in radio and the sharp rise in the dust spectrum or the 
dust peak near the rest wavelength of 100 /rm. Even when the dust peak in the FIR is 
not sampled by observations, the radio synchrotron measurements help set the vertical scale 
with respect to the dust spectrum, i.e. the SFR. 

Because SED measurements for most submm galaxies include only a few discrete points 
rather than a continuous frequency coverage, our best fit SED model search utilizes a min¬ 
imization with discrete sampling of the parameter space rather than a full cross-correlation 
technique. One advantage of this approach is that the upper limits in flux density can be in¬ 
corporated in a straightforward way by simply rejecting all trial SEDs that are incompatible 
with the upper limits. 


3.2. Trials on Galaxies of Known Redshifts 

To test the robustness of our SED template, we apply this photometric redshift tech¬ 
nique to several well studied submm galaxies, and these results are summarized in Table 2. 
Disregarding the formal uncertainties for the moment, the new photometric redshifts Zph are 
in excellent agreement with the spectroscopic redshifts Zgp for the all five submm galaxies 
with known redshifts. When compared to the old radio-to-submm spectral index estimates 
Zsi as shown in Figure 4, the improvement is seen mainly at high redshift {z > 2) where 
the effectiveness of the spectral index method diminishes due to the flattening of the a — z 
relation. Some improvement is generally expected for the new photometric method since 
more information is utilized. At the same time this comparison also highlights the efficiency 
of the spectral index technique which utilizes rather limited amount of information. 

The best fit models are shown in Figure 5 for the four cases for which the most SED 
data exist. For the two well studied SCUBA sources SMM J02399—0136 {z = 2.80, Ivison 
et al. 1998) and SMM J14011-I-0252 {z = 2.57, Ivison et al. 2000), we derive photometric 
redshifts of 2.83 & 2.73, respectively. The agreement with their spectroscopic redshifts is 
good in both cases, and the improvement over the a — z estimate is quite substantial for 



10 


SMM J02399—0136. The excellent agreement in SMM J02399—0136 is somewhat fortuitous 
since the 1.4 GHz measurement is well above the best £t SED model (thus the source of a 
large reduced but the best fit solution is limited correctly by the 3a upper limit at 8.7 
GHz. These two radio measurements are not consistent, and either the radio AGN is highly 
variable or the 1.4 GHz measurement may be affected by source confusion. The photometric 
redshifts {zph) for the submm galaxies SMM J02399—0134 {z = 1.06, Smail et ah 2000) 
and HRIO {z = 1.44, Gimatti et al. 1998; Dey et al. 1999) are also in good agreement with 
their spectroscopic redshifts within the uncertainty associated with the dust SED parameters 
(see § 4.4). The SED fit for HRIO shown in Eigure 5 suggests that the radio continuum is 
fainter than expected (i.e. fnth < 1), which drove the Zph to a slightly higher value. Only 
in GUDSS14.18 {z = 0.66, Lilly et al. 1999) does the photometric redshift differ from the 
spectroscopic redshift by an amount larger than the nominal uncertainty of the technique. 
Only one submm SED measurement exists for this galaxy, and this may be an important 
limiting factor. GUDSS14.18 also appears to be slightly underluminous in radio continuum. 


3.3. Photometric Redshifts for Submm Galaxies 

Since the great majority of the submm galaxies are too faint to yield optical red¬ 
shifts, a photometric technique may be the best way to infer their redshifts and evolu¬ 
tion. We apply our photometric redshift technique to seven submm galaxies whose spectro¬ 
scopic redshifts are unknown but are suspected of being at high redshifts: HDF850.1, Lock- 
man850.1, GUDSS14.1, FIRBAGK J1608+5418, SMM J00266+1708, SMM J09429+4658, 
SMM J14009-I-0252. The photometric redshifts and SFR we derive are summarized in Ta¬ 
ble 2, and the best fit models for six cases are shown in Figure 6. 

Perhaps the most interesting case is HDF850.1, which is the brightest SGUBA source in 
the Hubble Deep Field (Hughes et al. 1998; Downes et al. 1999a). Despite the extraordinarily 
deep HST imaging data available, this submm galaxy has eluded optical identification en¬ 
tirely. From the radio-to-submm flux density ratio, we have previously estimated its redshift 
to be 2 : > 2.6, but the flattening of the a — z relation at z > 2 meant that its redshift was not 
very well constrained (Garilli & Yun 1999). Our new photometric redshift analysis suggests 
that this galaxy has a high likelihood {Xn ~ 1) of having a redshift near ~ 4.1 with an 
intrinsic IR luminosity of about twice that of Arp 220. The uncertainty in the inferred red¬ 
shift is too large (a^ ~ 0.5) for a GO search using existing instruments, but a spectroscopic 
verification of HDF850.1 and other optically faint submm galaxies should become possible 
using future broadband instruments on the Green Bank Telescope or the Large Millimeter 
Telescope. 
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SMM J09429+4658 is another historically signihcant SCUBA galaxy, whose extremely 
dusty nature and red color (K = 19.4 mag, I — K >6) was clearly demonstrated for the first 
time by the combined analysis of the high angular resolution VLA continuum imaging and 
deep near-IR imaging (Ivison et ah 2000). Our new photometric redshift analysis suggests a 
moderately high redshift {zph ~ 3.9), but the SED model £t (Figure 6) is quite poor ~ 4) 
as our SED template cannot be fully reconciled with the three existing SED measurements. 
One possible explanation is that this submm galaxy is underluminous in radio continuum as 
in HRIO and CUDSS14.18, and its actual redshift is smaller. 

SMM J00266+1708 is another clear example of an optically faint {K = 22.5 mag) and 
red submm galaxy, whose extreme properties are clearly established by the combined analysis 
of the high angular resolution millimeter continuum observations at OVRO and deep near- 
IR imaging using the Keck telescope (Frayer et al. 2000). The derived photometric redshift 
of Zph = 3.50 makes it one of the highest redshift objects examined in this study, perhaps 
only the second after HDF850.1, and it has the second largest intrinsic FIR luminosity after 
SMM 02399-0136. 

CUDSS14.1 is the brightest SCUBA source discovered in the CFRS 14hr held by Eales et 
al. (2000). Its optical counterpart has been identified using high angular resolution imaging 
at the VLA and IRAM interferometer (Eales et al. 2000; Gear et al. 2000), but it has not 
yet yielded a spectroscopic redshift. All four significant detection points he along the model 
SED while all upper limits are consistent with the model. The derived photometric redshift 
of Zph = 2.06 is in good agreement with its optical photometric redshift oi z ^ 2 (Lilly et al. 
1999). Using the SED model of Dunne et al. (2000b) and optical to near-IR color. Gear et 
al. (2000) estimate its redshift to he between 2 and 4.5. 

Lockman850.1 is the brightest SCUBA source detected in the Lockman Hole ISOPHOT 
survey region with an extremely red {K ~ 20, 1 — K > 6.2) and extended optical counterpart 
(Lutz et al. 2001). Comparing the observed SED to that of Arp 220, Lutz et al. estimate 
its redshift to be around 3, which agrees well with our photometric redshift Zph = 2.72. 

SMM J14009-I-0252 is one of the submm sources identified in the Abell 1835 held by 
Ivison et al. (2000) with a faint optical counterpart {K ~ 21). We derive Zph = 1.30, but the 
SED £t shown in Figure 6 for the dust spectrum is poor {xi = 2.3). Ivison et al. rejected 
their own similarly modest redshift estimate because this galaxy is 3 times more luminous 
than HR 10 in the submm while it is more than 10 times fainter in K. Its 1.4 GHz hux 
density is rather large for a SCUBA source, indicating a possible presence of a radio AGN. 
The SED £t for the dust spectrum alone suggests 2 : ~ 3.5. 


FIRBACK J1608-I-5418 is one of the FIR sources identified by the FIRBACK deep 
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ISOPHOT survey (Dole et al. 2001). Incorporating its 350 /rm and 450 /rm detection at 
the CSO and the VLA detection at 1.4 GHz, Benford (1999) estimated its redshift to be 
around 2 : ~ 1.5. Our photometric redshift of Zph ~ 0.8 is substantially smaller. Its radio- 
submm flux density ratio suggests a presence of a radio AGN, and it is likely that higher 
than expected radio continuum flux density has led to a poor fit (y^ = 17). Excluding the 
radio measurements, an SED analysis for the submm/FIR dust peak alone suggests z ^ 0.7 
with about 4 times larger FIR luminosity than Arp 220. 


4. Discussions and Summary 

The above discussions of the SED modeling for the individual submm detected galaxies 
clearly indicate that systematic SED variation among the individual galaxies is the dominant 
factor over the statistical uncertainty associated with the SED measurements. Here we dis¬ 
cuss in detail the nature and the magnitude of the systematic effects associated with the dust 
and radio emission as well as the SED sampling and the intrinsic calibration uncertainties 
in the measurements. Their combined contribution to the photometric redshift uncertainty 
is quantihed in § 4.4. 

4.1. Effects of Dust Temperature and Emissivity 

Doppler shift of a dust spectrum to a higher redshift is exactly equivalent to lowering 
of dust temperature, and departure in from the template SED translates directly to a 
redshift error (see Blain 1999, for further discussions). In other words, if the dust temper¬ 
ature of a galaxy is actually higher than our template, then the resulting Zph would be an 
under-estimate. While inherently subject to an object-to-object variation, the magnitude of 
uncertainty due to the spread in is well understood quantitatively: ~ + z). If 

Td is higher or lower by 10 K, the resulting error in the photometric redshift is about 0.3 at 
2 : = 1, and it grows as (1 + z) reaching about 0.7 a.t z = 3. 

Dust emissivity f3 is as important as dust temperature Td because it determines the 
location of the low frequency edge of the thermal dust feature in the SEDs. There is signih- 
cant degeneracy between f3 and Td (e.g. Dunne & Eales 2001), and the decoupling the two 
becomes very difficult even when the dust spectrum is well sampled. In practice only a few 
measurements along the rising part of the dust SED are measured, and any dehciency in Td 
is at least partly offset by the compensating effect of (5. A trend supporting this effect is seen 
in Figure 2 as a broad trend of decreasing (3 with increasing Td- A better demonstration of 
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this effect is seen in Figure 1 where the submm SED points alone might be £t well by a dust 
spectrum with ~ 45 K and (3 ~ 1.5. The observed trend of increasing and deceasing 
(3 with increasing SFR is also clearly seen in the theoretical modeling of dust heating and 
emission using a self-consistent 3-D radiative transfer code by Misselt et al. (2001). Therefore 
the error estimate based on the Td alone might serve only as a reasonable upper bound. 


4.2. Effects of the Radio Spectrum 

Another reason why the photometric redshift error based on the scatter in Td alone 
might be an over-estimate is that our SED fitting scheme relies as much on the spectral 
trough between the radio continuum and thermal dust emission as the dust spectrum itself. 
If this spectral trough is considered as the primary SED feature from which the photometric 
redshift information is derived, then it is easy to see that the variation in the dust spectrum 
is strongly modulated by the radio spectrum which is not affected. 

The radio spectrum has its own uncertainty in its overall scaling as shown by the dis¬ 
tribution of fnth (see Figure 3), and its contribution to the photometric redshift error has 
already been alluded to in § 3.2. Presence of a powerful radio AGN was a severe limiting 
factor for our earlier redshift estimation based on the radio-to-submm spectral index (Carilli 
& Yun 1999), and it is not surprising that radio AGNs pose the biggest obstacle for the SED 
htting photometric technique as well. The photometric redshifts derived for an ensemble 
of QSOs that are associated with dusty hosts are given at the bottom of Table 2 (also see 
Fig. 4). Although there are exceptions such as LBQS 1230-1-1627 where Zph ~ Zgp, the new 
photometric method does not fare much better than the spectral index method - the median 

Xl ~ 7 . 

In practice, radio AGNs are not likely to limit the photometric redshift technique utiliz¬ 
ing a starburst SED template because only a small minority of the known submm galaxies 
show evidence of hosting a powerful AGN. Gomparisons of deep X-ray and submm surveys 
have shown very little overlap between the detected sources (e.g. Fabian et al. 2000; Horn- 
schemeier et al. 2000; Barger et al. 2001). At least three out of five submm galaxies discussed 
in § 3.2 are known to be an AGN host, but the presence of an AGN also makes their spectro¬ 
scopic redshift measurement possible and the use of a photometric technique unnecessary. 
The total fraction of infrared selected galaxies whose radio-to-FIR SED shows a clear sign of 
energetic AGN is only a few percent in the local universe (Yun et al. 2001), and this fraction 
appears to remain roughly the same at higher redshifts. 
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4.3. SED Sampling and Calibration 

Another important source of uncertainty is poor sampling of the SED. Sampling only 
the radio part of the SED offers virtually no redshift information, and the same is true if only 
the submm part of the SED is measured. On the other hand, because we are htting the SED 
features changing over logarithmic scales, only a few well placed SED data points are needed 
to derived the redshifts (Carilli & Yun 1999). More than one measurement on both side of 
the SED trough is highly desirable for a successful photometric redshift determination. 

A related problem is the calibration and the relative weighting of individual SED mea¬ 
surements. Measurement uncertainties found in the literature often do not fully account for 
the overall calibration accuracy as demonstrated by the magnitude of the scatter and the 
sizes of the error bars plotted for the SED measurements of Arp 220 in Figure 1. For exam¬ 
ple, the IRAS 60 /rm and 100 /rm measurements found in the Point Source Catalog (PSC) 
and the Faint Source Catalog (ESC) are known to have systemic calibration differences of 
order 10% while the flux densities reported in the ESC often include uncertainties of only 
a few percent. For deriving the SED template from the 23 infrared selected starbursts, we 
adopted flux densities in the ESC but increased measurement uncertainties to 10% (unless 
the reported uncertainty was larger). This re-calibration of the measurement uncertainty is 
more than cosmetic since the relative weight of the data points are directly reflected in the 
minimization. We found this re-weighting to be critically important in utilizing the submm 
measurements. Similarly, all radio continuum measurements are assigned at least 10% un¬ 
certainty in order to account for the overall uncertainty in the flux calibration. No effort 
was made to re-calibrate any of the submm data points since most measurements generally 
carry relatively large fractional uncertainties, but some are clearly under-estimates based 
on the SED plots such as shown in Figure 1. The same problems also plague the reported 
SED measurements of many submm galaxies, and they will inevitably impact the accuracy 
of the photometric redshifts derived from these data. One way to improve the situation in 
the future would be employing a frequency selective bolometer (e.g. Kowitt et ah 1996; 
Meyer et al. 2001) that can make simultaneous measurements of several submm bands with 
accurate relative calibration between the measurement bands. 


4.4. Estimate of Uncertainty in the Radio-to-FIR SED Technique 

The uncertainty of az ~ 0.3(1 -|- z) derived from the scatter in dust temperature is 
probably the upper bound to the error one may expect from the photometric technique 
using the dusty starburst SED template. The compensating effects of the dust emissivity 
and the radio continuum are more difficult to quantify as are the uncertainties contributed 
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by the scatter in fnth- One way to estimate the collective uncertainty of the SED photometric 
technique is to apply this technique to the same 23 dusty starbursts from which the SED 
template was derived. This is not entirely circuitous since our SED template is based only 
on the average properties of these galaxies and has no knowledge of dispersions. At the 
least, we may be able to conhrm the impact of the 10 K scatter in the dust temperature on 
the dispersion in Zph if dust temperature variation dominates the uncertainty in determining 
photometric redshifts. 

The resulting “photometric redshift” Zph for the 23 dusty starbursts are listed in the 
last column of Table 1. There are several galaxies with negative Zph since the impact of 
higher dust temperature for a 2: = 0 galaxy with > 58 K would result in Zph < 0. The 
minimization program is modihed to search a redshift range of — 1 > Zph > +1 with modihed 
scaling along the flux density axis in order to remove the non-physical impact of negative 
redshift on Dl- A handful of galaxies with negative redshift are known (e.g. M81), but 
negative redshift is generally considered non-physical. Here, the only physically meaningful 
interpretation of a negative Zph we derive is the measure of the magnitude of departure in 
the radio and dust property from the template SED, just the same way any uncertainty in 
Zph should be interpreted at any redshift. 

A histogram of the resulting photometric redshifts shown in Figure 7 suggests that the 
median “redshift” is about -1-0.05, suggesting the template SED may be slightly biased, but 
this offset is probably not very signihcant given the range of “redshifts” derived. “Photomet¬ 
ric redshifts” as large as -1-0.3 and as small as —0.3 are found as expected from the scatter 
in Td- The two extreme negative ‘Tph” object IRAS 08572-1-3915 and Mrk 231 have Td of 74 
K and 72 K, respectively, following the general trend expected of the dust temperature vari¬ 
ation. On the other hand, two other dusty starbursts with characteristic Td > 70 K, IRAS 
15250-1-3609 and IRAS 05189—2524, have ‘Tph” of -fO.Ol and -1-0.05, clearly demonstrating 
that there are other compensating effects and that the scatter in Td alone does not dictate 
the overall photometric redshift uncertainty. 

Regardless of the underlying causes, the “photometric redshifts” for 2/3 of all galaxies 
he within Az < 0.10, and we estimate the collective uncertainty of the photometric redshift 
technique, including the variations in the radio and dust properties as well as the uncer¬ 
tainties in the minimization process, is about 0.1 as long as this technique is applied to 
luminous dusty starburst galaxies only. Allowing for the (1 -|- ;2) frequency folding of the 
Doppler effect, we estimate an overall uncertainty of ~ 0.1(1 -|- z) with an upper bound 
in redshift uncertainty of about 0.3(1 + 2:). 

In either case, this photometric redshift technique utilizing the radio-to-FIR dusty star- 
burst SED represents a signihcant step forward, particularly at high redshifts, when com- 
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pared with existing methods. The full potential of this method will be realized when the 
sources identihed by several large deep, multi-frequency surveys planned in the immediate 
future (e.g. SIRTF Legacy Surveys) are analyzed together to reveal an accurate redshift 
distribution of luminous dusty galaxies at high redshift. 


The National Radio Astronomy Observatory is a facility of the National Science Foun¬ 
dation operated under cooperative agreement by Associated Universities, Inc. Some of the 
data presented here are obtained from the NASA/IPAC Extragalactic Database (NED), 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology, un¬ 
der contract with the National Aeronautical and Space Administration. 
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Fig. 1.— Radio-to-IR SED for Arp 220 is shown along with the model starburst SEDs with 
(3 = 1.35 and = 43 K, 58 K, and 73 K are shown. The radio-FIR normalization term fnth 
is set to 1.0, which is the Galactic normalization by Condon (1992). 
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Fig. 2.— A distribution of best fit (3 values versus Td for the 23 luminous IR starb 
galaxies. The median values for (3 and Td are 1.35 & 57 K while the mean values 
1.32 ± 0.17 and 58 ± 9 K, respectively. 
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Fig. 4.— (a) A comparison plot of spectroscopic redshifts (zopt) versus redshift estimates from 
the radio-to-submm spectral index {zsi). (b) A comparison plot of spectroscopic redshifts 
{zopt) versus radio-to-FIR SED photometric redshifts (zph). Solid squares represent submm 
galaxies while the empty circles are submm detected optical QSOs. 
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Fig. 5.— Results of the SED fits for four sub mm galaxies with known spectroscopic redshifts. 
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Fig. 6.— Results of the SED fits for six sub mm galaxies with unknown redshifts are shown, 
including HDF850.1. 











Table 1. Spectral Energy Distribution Model Fits 


Name 

L(FIR) 

(10'°Le) 

SFR^ 
{Mq yr-^) 

Td 

(K) 

/9" 

fnth'' 

''Zph” 

Zw049.057 

10 

18 

54 

1.30 

0.6 

+0.06 

NGC 5256 

11 

19 

51 

1.35 

4.0 

+0.07 

MCG+00-29-023 

11 

19 

49 

1.15 

1.3 

+0.25 

NGC 3110 

12 

20 

46 

1.35 

1.2 

+0.19 

Mrk 331 

16 

27 

57 

1.15 

0.8 

+0.06 

NGC 1614 

21 

37 

61 

1.40 

1.0 

-0.07 

UGC 2369 

21 

37 

52 

1.30 

1.2 

+0.16 

NGC 2623 

21 

37 

57 

1.45 

1.1 

-0.07 

Arp 193 

25 

42 

47 

1.70 

1.5 

+0.01 

Arp 148 

25 

42 

50 

1.20 

0.9 

+0.22 

IRAS 10173+0828 

33 

57 

69 

1.10 

0.6 

+0.02 

NGC 6240 

35 

61 

49 

1.45 

3.9 

-0.07 

Mrk 848 

35 

61 

62 

1.20 

1.6 

+0.06 

IRAS 15250+3609 

48 

82 

74 

1.05 

1.1 

+0.01 

UGC 5101 

54 

93 

47 

1.40 

2.5 

+0.14 

IRAS 08572+3915 

55 

93 

74 

1.60 

0.3 

-0.22 

IRAS 05189-2524 

57 

98 

70 

1.05 

0.6 

+0.05 

IRAS 10565+24 

60 

104 

55 

1.55 

0.4 

+0.06 

Mrk 273 

73 

127 

61 

1.40 

1.7 

-0.10 

Arp 220 

95 

163 

59 

1.15 

1.0 

+0.05 

IRAS 12112+0305 

106 

183 

63 

1.25 

0.8 

+0.02 

IRAS 14348-1447 

109 

187 

63 

1.50 

1.0 

+0.09 

Mrk 231 

137 

236 

72 

1.35 

2.0 

-0.30 


‘^SFR = LFiji/{5.8 X IO^Lq) Mq yr ^ (Kennicutt 1998) 

'’Dust emissivity. 

'’Multiplicative correction factor for non-thermal radio continuum emis¬ 
sion. The fnth = 1-0 corresponds to the Galactic normalization adopted by 
Condon (1992). 
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Table 2. Redshift Estimates for Submm Sources 


Name 

^sp 

zsi 

^ph 

SFR^ 

Ntot {Naet)'° 

2 c 
An 

Refs. 

Known Redshifts: 








CUDSS14.18 

0.66 

1 11^+0.53 
-L--L^_0.45 

1.11 + 0.21 

200 

5 (3) 

0.02 

1 

SMM J02399-0134 

1.06 

1 1 'T+O.SS 
-0.45 

0.93 + 0.19 

525 

5 (3) 

0.5 

2 

HRIO 

1.44 

J--4y_0.53 

1.75 + 0.28 

380 

7 (5) 

0.7 

3,4,5 

SMM J14011+0252 

2.57 

9 ko+1-24 
^•O0_o g2 

2.73 + 0.37 

850 

5 (4) 

1.0 

2,6 

SMM J02399-0136 

2.80 

J-.OO_0 53 

2.83 + 0.38 

1400 

9 (5) 

28 

2,7 

No Known Redshifts: 








FIRBACK J1608+5418 

- 

U.UO_Q Q5 

0.84 + 0.18 

690 

7 (5) 

17 

12,13,14 

SMM J14009+0252 

- 

1 q't+0.59 
-*-•^^-0.54 

1.30 + 0.23 

980 

5 (4) 

2.3 

2,6 

CUDSS14.1 

- 

2 01 

2.06 + 0.31 

520 

7 (4) 

0.5 

10,11 

Lockman850.1 

SMM J00266+1708 

— 

Z,.yO_Q gg 

o.^y_i 23 

2.72 + 0.37 

3.50 + 0.45 

570 

1200 

8 (5) 

6 (4) 

0.2 

0.2 

15 

2 

SMM J09429+4658 

- 

> 3.6 

3.86 + 0.49 

890 

4(4) 

4.0 

2,6 

HDF850.1 

- 

> 2.6 

4.11 + 0.51 

380 

8 (4) 

0.8 

8,9 

QSOs: 








FSC 10214+4724 

2.28 


1.51 + 0.25 

2600 

9 (9) 

7.6 

16,17 

H 1413+117 

2.56 

0.88t°;" 

0.31 + 0.13 

850 

9 (9) 

28 

17,18 

LBQS 1230+1627 

2.74 

2 19+^-^® 
^•J^W-0.75 

2.66 + 0.37 

1300 

4 (4) 

0.5 

19,20,21 

SMM J04135+1027 

2.83 

> 0.9 

1.83 + 0.28 

1560 

4 (3) 

0.8 

22 

APM 08279+5255 

3.87 

9 1 o + l-lO 

^•lo_o.7i 

1.74 + 0.28 

4900 

7 (7) 

9.4 

23,24,25 

BR 1335-0417 

4.40 

1 oi +0.95 

l-ol_0.62 

2.84 + 0.38 

1200 

5 (5) 

7.0 

17,19,20,21,26 

BR 0952-0115 

4.43 

1 qi::+0.63 

1.00_Q 4g 

0.96 + 0.20 

460 

5 (4) 

5.6 

19,20,21,26 

BR 1202-0725 

4.70 

O 1 0 + 1-81 

'^•-‘-'^-1.14 

3.20 + 0.42 

2300 

11 (9) 

3.1 

17,19,20,21,26 


^SFR is in Mq yr 

^Ntot and N^et are total number of SED data used for the fit and the number of detections (i.e. Ntot = 

^det ^limits)' 

’'Xn = X^l^det- 

References. — (1) Eales et al. (2001); (2) Smail et al. (2000); (3) Cimatti et al. (1998); (4) Dey et al. 
(1999); (5) Andreani et al. (2000); (6) Ivison et al. (2000); (7) Ivison et al. (1998); (8) Hughes et al. (1998); 
(9) Downes et al. (1999a); (10) Eales et al. (2000); (11) Gear et al. (2000); (12) Benford (1999); (13) Scott 
et al. (2000); (14) Dole et al. (2001); (15) Lutz et al. (2001); (16) Rowan-Robinson et al. (1993); (17) 
Benford et al. (1999); (18) Barvainis et al. (1992); (19) Omont et al. (1996); (20) Guilloteau et al. (1999); 
(21) Yun et al. (2000); (22) Knudsen et al. (2000); (23) Irwin et al. (1998); (24) Lewis et al. (1998); (25) 
Downes et al. (1999b); (26) McMahon et al. (1999) 



